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We analyze the magnitude of the specific heat jump AC at the superconducting transition tem- 
perature T c in the situation when superconductivity develops in the pre-existing antiferromagnetic 
phase. We show that AC/T C differs from the BCS value and is peaked at the tri-critical point where 
this coexistence phase first emerges. Deeper in the magnetic phase, the onset of coexistence, T c , 
drops and AC/T C decreases, roughly as AC/T C oc T c 2 at intermediate T c and exponentially at the 
lowest T c , in agreement with the observed behavior of AC/T C in iron-based superconductors. 

PACS numbers: 74.25.Bt, 74.25. Dw,74.62.-c 



Introduction. The magnitude and the doping depen- 
dence of the specific heat jump at the superconducting 
transition temperature T c is one of unexplained phenom- 
ena in novel iron-based superconductors (FcSCs) [l|. In 
BCS theory AC/T C ~ 1.437, where 7 = ir 2 N F /3 is the 
Sommcrfcld coefficient, and Np is the total quasiparti- 
clc density of states (DoS) at the Fermi surface (FS) in 
the normal state. Although the behavior of FcSCs is in 
many respects consistent with BCS theory, the exper- 
imental values of AC /T c vary widely between different 
compounds, ranging between 1 mJ /(mol ■ K 2 ) in under- 
doped Ba(Fe 1 - x Ni x ) 2 As 2 and 100 mJ/(mol ■ K 2 ) 
in optimally hole-doped Ba\- x K x Fe2As2 [3j. Such 
huge variations may be partly due to differences in 7, 
which were indeed reported to be larger in hole-doped 
FeSCs 0, 0). Yet, even for a given material, e.g., 
Ba(Fei- x Ni x )2As2 or Ba{Fe\- x Co x )2As2 the magni- 
tude of AC/T C is peaked near optimal doping x op t and 
rapidly decreases, approximately as AC /T c oc T 2 , at 
smaller and larger dopings 0- 

This rapid and non-monotonic variation of AC /T c over 
a relatively small range of 0.03 < x < 0.12 is unlikely to 
be attributed to change in 7 and has to be explained by 
other effects. According to one proposal @, AC/T C oc T 2 
is caused by interband scattering off non-magnetic impu- 
rities, which is pair-breaking for pairing. However, 
AC/T C oc T? only holds in the gapless superconductivity 
regime, 0, [3| when T c is already reduced by more than 
factor of five as compared to T c for a clean case. Be- 
fore that, the reduction of AC/T C by impurities is rather 
mild. Experimentally, the reduction is strong immedi- 
ately away from optimal doping and, moreover, occurs 
on both sides from the optimal doping. This reduction 
is difficult to explain by impurity scattering. Large value 
of ACJT C at x opt may also be due to strong coupling ef- 
fects |J| . This is certainly a possibility, but we note that 
AC/T C is a non-monotonic function of the coupling Q 
and at strong coupling is actually smaller than the BCS 
value. 

We propose a different explanation. We argue that 



the origin of strong doping dependence of AC/T C is the 
coexistence of spin-density-wave (SDW) magnetism and 
s ± superconductivity (SC). 0-Q In Ba(Fe\- x Ni x ) 2 As2, 
Ba{Fe\- x Co x )2As2 and, possibly, in other FeSCs, op- 
timal doping x op t nearly coincides with the end point 
of the coexistence region (tri-critical point). Q We an- 
alyzed the behavior of AC/T C near x op t within a mean- 
field BCS-like theory and found that AC/T C is by itself 
discontinuous and jumps by a finite amount when the 
system enters the coexistence region, see Figs. [1] and [2] 
The magnitude of the jump depends on microscopic pa- 
rameters of the system and for a wide range of parame- 
ters, AC/T C exceeds the BCS value. Beyond a mean-field 
treatment, paramagnetic fluctuations transform the dis- 
continuity in AC/T C at x op t into a maximum, such that 
AC /T c smoothly decreases on both sides of optimal dop- 
ing x pt , as illustrated in Fig. [TJ 

We also examined the behavior of AC/T C along the 
entire transition into the coexistence state. We found 
that AC/T C decreases together with T c , and eventually 
becomes smaller that its BCS value. Deeper in the SDW 
region, AC/T C oc T 2 and vanishes exponentially at the 
lowest T c , see Fig. [2] The explanation of this behavior 
goes beyond a standard paradigm that T c and AC de- 
crease because FS available for SC is modified by SDW. If 
that was the only effect, then the DoS would not change 
significantly and AC/T C would only weakly depend on 
T c . We found, however, that the SC transition line nec- 
essarily enters the region in which SDW order gaps out 
the whole FS. In this region, SC appears because the sys- 
tem finds that it is energetically advantageous to decrease 
the magnitude of the SDW order parameter and develop 
a non-zero SC order parameter, even in the absence of 
the FS. Because all low-energy states in this region are 
gapped at T c , the thermodynamic characteristics, includ- 
ing AC/T C , are exponentially small. We found that the 
precursors of this behavior develop at higher T c , when 
the reconstructed FS is still present. As a result, AC/T C 
decreases as the SC temperature, T c , drops. 

The behavior of AC/T C outside the coexistence region 
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FIG. 1: (Color online) (a) The specific heat C(T) and SC and 
SDW order parameters A and M as functions of T. We con- 
sider the jump of C(T) at the onset of SC. (b) The behavior 
of AC /{^T c ) as a function of <5o, which scales with doping. In 
a mean-field theory, AC/T C is discontinuous at the end point 
of the coexistence state (P) and jumps back to BCS value 
AC/T C ~ 1.437 a t larger dopings (dashed horizontal line). 
Beyond mean-field, paramagnetic fluctuations smear the dis- 
continuity of AC/T C and transform it into a maximum, as 
schematically shown by the solid line. 



is likely to be a combination of several effects. When 
paramagnetic fluctuations weaken, AC /T c reduces to its 
BCS value. Further decrease of AC/T C is partly due to 
impurities, [ll, [l3[ and partly due to shrinking of the hole 
FSs and to the fact that at larger x the gap along electron 
FSs becomes more anisotropic. 

The method. To obtain AC, we expand the free en- 
ergy in powers of the SC order parameter A to order A 4 . 
When the SC transition occurs from a pre-existing SDW 
state, the expansion reads 
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FIG. 2: (Color online) Top: The phase diagram in T-So plane 
for T mfi /T cfi = 2 and several 6 2 /(2itT c , ) = 0.4, 0.28, 0.26, 
corresponding to wide, medium, and narrow doping ranges of 
the coexistence phase. SDW, SC and the SDW+SC phases 
meet at the tri-critical point P (in this case also meeting nor- 
mal (N) state). Bottom: The behavior of AC/-yT c vs T c /T c , 
in the coexistence region for these 82- The arrows indicate 
T c , below which the whole FS is gapped by SDW. As T c is 
lowered through this value, AC/T C decreases, as T c 2 at inter- 
mediate T c and exponentially at lower T. Near the tri-critical 
point, AC /Tc may well exceed the BCS value 1.437- 
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k x and k y are deviations from (it, it). The same model 
has been considered in Refs. |8l4lu| on the coexistence of 
SDW and SC orders. At perfect nesting, £ e = —£h, while 
a A (M , T)A 2 + ri(M , T)A 4 , (1) for a non-perfect nesting £ e = -^+26^ near the FS, and 

Stp = Sq + 62 cos 2ip captures the difference in the chemical 



where J-q = J-(0, Mo) is the free energy of a pre-existing 
SDW state, M = M (T) is the SDW order parame- 
ter which minimizes J-(0,M), and 77 includes the feed- 
back of the finite SC order parameter on the SDW state, 
A/ 2 = A/ 2 - C(A 2 ), see Fig. CfTa). The T c is given by 
a/±(M (T c ), T c ) = and the specific heat jump is 



<9cva dMy 
dM 2 ~dF' 



AC 37 / dctA \ 2 dctA daA 

. . (2) 

As rj — > 0, A appears discontinuously and the transition 
between SDW and SC becomes first order. 

To obtain actual expressions for cva and rj, we need to 
specify the band structure of a material. Since our goal 
is to demonstrate the discontinuity of AC /T c at x op t and 
the reduction of AC/T C along the coexistence onset, we 
adopt a simplified 2D two-band model with the hole-like 
band near the center of the Brillouin zone (BZ), with 
£,h = /i/i-fc 2 /2m/„ and electron-like band near the corner 
of the BZ, with £ e = — /i e + k 2 /2m x + k 2 /2m y , here 



potentials and in electron m x _ y and hole masses via 
(5o , and anisotropy in m x and m y via 82- Without loss 
of generality, we assume that 80 changes with doping, 
but the ellipticity parameter 62 remains intact. We use 
a model with four-fermion interactions in SDW and SC 
channels 0, [ll], 12 1, decompose them using SDW and 
SC order parameters M and A, and express couplings 
in terms of transition temperatures T Ci o to the SC state 
in the absence of SDW and T m .o to the perfectly nested 
SDW state, <5o,2 = 0) in the absence of SC. Note that the 
actual T m decreases, even in the absence of SC, when 80 
and 82 increase, while T c $ is the transition temperature 
at the tri-critical point indicated by P in Fig. [2] 
The free energy for such a model has the form [j| 
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where E n = ^e 2 + A 2 , e n = nT{2n + 1) arc the Mat- 
subara frequencies (n = 0, ±1, ±2, . . . ), and (...) denotes 
averaging over tp along FSs. For this functional we find 
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The derivatives are taken at A = and M = Mo with 
Mq defined by 
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In the absence of SDW, M = 0, da A /dT = da A /dT, 
T] = A(M = 0) = 7C(3)/(32tt 2 T 2 ), and we reproduce 
the BCS result AC/T C = 1.43j. To obtain AC/T C inside 
SDW phase we solve Eq. fl5j for M 2 (T), insert the result 
into Eqs. evaluate daA/dT and ?y and finally use 
Eq. ([2]). AC/T C depends on three input parameters So, S2, 
and T rn ,o/r Ci o, and generally differs from the BCS value. 

Results. We present AC/T C as function of So for fixed 
S2 and T m fi/T C fi in Fig.QJb). It grows from zero value at 
the low-temperature onset of the coexistence phase and 
reaches its maximum at the tri-critical point, where SDW 
order disappears. At this So, T c reaches T Cj o and AC/T C 
jumps to the BCS value. 

Plotted as a function of T c /T cfi in Fig. [2 AC/T C shows 
exponential behavior at small T c and approximate T 2 be- 
havior at intermediate T c /T c ,o < 0.5. The magnitude 
of AC/T C at T Cj o increases when the width of the coex- 
istence region shrinks. This can be easily understood, 
since shrinking of the coexistence region brings the sys- 
tem closer to a first order transition between SDW and 
SC at which the entropy itself becomes discontinuous at 
T c , and AC/T C diverges. In the opposite limit, when the 
width of the coexistence range is the largest, the magni- 
tude of AC/T C is smaller and can even be below the BCS 
value. 

The behavior of AC/T C at small T c and at T c w T c . 
can be understood analytically. We first focus on the low 
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FIG. 3: (Color online) Left: the coexistence region (un- 
shaded) in the 82-80 plane. Each point corresponds to a 
particular ratio T m fi/T c fi and this ratio increases monoton- 
ically as So grows at fixed 82. Right: The value of AC/-yT c 
at the end point of the coexistence region, at T c = T c .o — 
for S 2 /2irT c , = 0.4,0.28,0.2 Thin solid line is the BCS value 
AC /T c = I.437. Over some range of parameters, AC/"/T c at 
Tc.o — significantly exceeds the BCS value. Diamonds rep- 
resent the values AC /^T c at T c — T c .q — for the curves for 
S 2 /27vT cfi = 0.4 and 0.28 in Fig. [2] 



T c limit, T c <C T Ct o- We argue that at low T c a conven- 
tional reasoning that the mixed SDW/SC state exists be- 
cause FS reconstruction associated with SDW order still 
leaves pieces of the FS available for SC does not work 
because the SDW state immediately above T c is already 
fully gapped. Indeed, to a logarithmical accuracy, the 
condition aA(Mo,T c ) = becomes 



27rT c — ( 1 - Re 

e n >0 £n \ 




= ln 



T, 
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which is only satisfied if Af > max{5 v } = S a + <5 2 , i.e. 
the SDW state is fully gapped [§]. 

The coexistence state emerges from the fully gapped 
SDW state because, when 77 > and qa > 0, it is en- 
ergetically advantageous to gradually reduce the magni- 
tude of the SDW order M below Mo, thus change sign 
of «a and create non-zero SC order A. The contribution 
to the free energy from SC ordering comes from the re- 
arrangement of quasiparticle states above the gap, and 
the number of such quasiparticles exponentially decreases 
with lowering T c . As a result, the magnitude of the spe- 
cific heat discontinuity at T c also becomes exponentially 
small, AC/T C oc exp[— (M — So — 62) /T c ], consistent with 
Fig. [2] Wc emphasize that as long as SDW-coexistence 
transition is of second order, r\ > 0, the exponential be- 
havior at low T c always holds. The behavior of AC /T c 
at larger T is less universal, but the reduction of AC /T c 
always begins even before the condition Mo > So + 82 is 
reached, see Fig. [2] For smaller r\, AC/T C decreases with 
T c over the whole range T c < T c o- 

We next consider AC/T C near the end point of the co- 
existence regime, when T m — > T Cj o + 0, T c — > T c ,o — 
and Ma is small. In this case, we expand a a, i], 
and Jo in terms of M 2 , use <9a A (0, T)/dT = 1/2T, 
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da A /dM 2 = 2C and dM 2 /dT = -(da m /dT)/2B and 
express Eq. © as 

AC = 3 7 (1 da m C \ 2 
T c 27r 2 77 \2T dT 5 / T=Tc ' [ ' 

where 770 = ^0 — Cq/Bq and the coefficients Aq, Bq and 
Co are given by Eqs. (gj with M = 0. For dM^/dT 
we obtain from Eq. © dM$/dT = -(da m /dT)/2B , 
Eq. with 

1 2 ^ / \ rs^ 

2T - 2-, o \( £ 2+«^/- 1 j 

In the absence of SDW order, AC/T C does not contain 
(C Q /B )da m /dT, r, Q -> A , and AC/T C = I.437. Once 
M is small but finite, AC/T C is determined by the inter- 
play between (C /B )da m /dT and Cq/B in 779- Both 
terms originate from the fact that SDW order is sup- 
pressed by SC order, M 2 = M§ - (C/B)A 2 . As a result, 
AC/T C is different from the BCS value already for in- 
finitesimally small SDW order and exhibits discontinuity 
upon crossing the SDW border. In general, the change 
in 770 is more important because C]: — AqBq at perfect 
nesting, 5 ,2 = 0, and 170 =0 [a ls[- Reduced value of 
77 leads to generally larger AC/T C at T c = T c , — than 
AC/T C = 1.437 f° r transition from normal state to a pure 
SC state. As the coexistence region narrows, 770 — > and 
AC/T C grows. Figure [5] shows AC/T C at T c = T cfi - for 
770 > 0. In a wide range of parameters, AC /T c > I.437. 

Beyond mean-field. In a mean- field description, AC/T C 
is discontinuous at the tri-critical point P, Fig. [2J with 
T m = T c . The free energy F and the specific heat 
jump depend on the finite square of the SDW order, 
M§ oc (T m - T). Although above T m the average 
Mo = 0, one expects to replace M5 2 by the finite sec- 
ond moment of SDW order due to Gaussian fluctua- 
tions, (A/5 2 } fi uct oc (T — T m ). These fluctuations enhance 
AC /T c and transform the discontinuity in AC /T c into a 
maximum, as shown in Fig. [T] As a result, AC/T C drops 
for deviations from T C: o = T m both into the coexistence 
phase and away from the SDW region. Still, the decrease 
of AC/T C should be more rapid within the SDW-ordered 
phase. An enhancement of AC/T C by paramagnetic fluc- 
tuations was earlier obtained in Ref. [14|. 

Conclusions. We demonstrated that the specific heat 
jump AC/T C across transition from SDW to the co- 
existence phase significantly deviates from the BCS 
value. AC/T C has its maximum for doping at the high- 
temperature end of the coexistence phase and decreases 
for doping deviations in both directions. In the co- 
existence phase, AC/T C eventually becomes exponen- 
tially small as all low-energy states at low T c are al- 
ready gapped out by the SDW order. At intermediate 
T c ^ T C fi, AC/T C scales approximately as T 2 . 

This behavior is quite consistent with the observed 
doping evolution of AC/T C in Ba(Fe\- x Ni x )2As2 and 



Ba{Fe\- x Co x )2As2 0- in these materials AC/T C is 
peaked at the tri-critical point, which coincides with the 
optimal doping x opt , and decreases for deviations from 
x op t in both directions, faster into the coexistence region. 
The evolution of AC/T C there approximately follows T 2 
behavior. 

As doping increases and paramagnetic fluctuations dis- 
appear, AC/T C reduces to the BCS value. Further reduc- 
tion of AC/T C at larger dopings x > x opt is, most likely, 
a combination of several effects: (1) the enhancement of 
a non-magnetic interband impurity scattering (2) 
stronger anisotropy of the gap on the electron FSs that 
increases 77; (3) the reduction of 7 due to shrinking of the 
hole FSs. 

How strongly the value of AC/T C at x op t exceeds the 
BCS result is difficult to gauge because 7 has to be ex- 
tracted from the normal state C(T) for which 7T contri- 
bution is only a small portion of the total specific heat. 
In Ba{Fe 1 - x Co x )2As2, AC/T C ~ 26 mJ/(mol ■ K 2 ) at 
Xopt, an d 7 — 20 mJ/ (mol-K 2 ) In this case, the max- 
imum of AC/T C is not far from the BCS result, and the 
reduction of AC/T C into the coexistence state is mainly 
due to the gapped low-energy states by SDW order at 
low T c . At the same time, a very strong AC/T C over 
a 100 mJ/(mol ■ K ) near x op t in Ba\- x K x Ft2As2 is 
well above the BCS value 0, even if 7 is as large as re- 
ported [1,0] 50 — 60 mJ/(mol ■ K 2 ). This enhancement 
can well be due the peak in AC/T C at the onset of the 
coexistence regime or strong SDW fluctuations. 
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